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Abstract 

We review the development of renormalizable SO(IO) Susy GUTs based on the 
210 © ©126 © 126 © 10 Higgs system. These GUTs are minimal by parameter 
counting. Using the 50(10) SU{4) x SU{2)l x SU{2)r label decomposition 
developed by us we calculated the complete GUT scale spectra and couplings and 
the threshold effects therefrom. The corrections to ac, Sin^Ow and Mx are sensi- 
tive functions of the single parameter ^ that controls symmetry breaking and slow 
functions of the other parameters. Scans of the parameter space to identify regions 
compatible with gauge unification are shown. 

The tight connection between the phenomenology of neutrino oscillations and 
exotic (AS / 0) processes predicted by Susy SO(IO) GUTs is discussed in the con- 
text of the recent successful fits of all available fermion mass/mixing data using the 
10 © 126 Higgs representations and Type I/II seesaw mechanisms for neutrino mass. 
We emphasize that the true output of these calculations should he regarded as the 
unitary matrices that specify the orientation of the embedding of the MSSM within 
the MSGUT. The AB ^ 0, d=5 effective Lagrangian is written down using these em- 
bedding matrices. We show how Fermion Mass fitting constraints can be combined 
with GUT spectra to falsify/constrain the MSGUT and its near relatives. An initial 
survey indicates that Type I Seesaw neutrino masses dominate Type II and both 
are too small in the perturbative MSGUT even when the mixing and mass squared 
splitting ratios are as per data. This motivates a detailed study of the MSGUT con- 
straints using the outputs of the fitting of fermion data, as well as consideration of 
modifications/extensions of the MSGUT scenario. 
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1 Introduction 



The Supersymmetric SO(IO) GUT based on the 126 © 126 © 210 Higgs multiplets was 
first introduced over 20 years ago in^|2]. At that time the electroweak gauge parameters 
and fermion mass data were quite incomplete. The next significant development was the 
realization |2j that R parity Rp = (_ )3(B-i)+25 formed a part of the gauge symmetry in LR 
symmetric models with gauged f/(l)B-L[4 and was hence protected by it. In particular it 
would remain unbroken as long as no field with odd B — L received a vev. SO(IO) is now 
the most favoured GUT gauge group because of the natural way in which it accommodates 
complete fermion families together with the superheavy right handed neutrinos required by 
the seesaw mechanism In 1992, with LEP data in hand and CKM parameters largely 
known Babu and Mohapatra|ni took up the task of fitting Fermion masses and mixing using 
the 10 © 126 FM Higgs system introduced injUEj and the Type I seesaw mechanism. They 
proposed that once the charged fermion masses in the SM had been fitted the 10 © 126 
matter fermion Yukawa couplings of the GUT would be completely determined and hence 
the model would predict the allowed neutrino masses. Although, based on the then current 
notions about neutrino masses/mixings, they concluded that their proposal failed, it would 
inspire much future work [71 |H1 El UHl HI]- 

In the mid-1990's the Supersymmetric LR model based on a parity odd singlet was 
found to havep^ a charge breaking vacuum and this was construed as evidence of a low 
scale linked breaking of SU{2)r, B — L and R-parity[T3] driven by the soft supersymmetry 
breaking terms. However an alternative analysis based on allowing this scale to lie anywhere 
in the range from Ms and Mx allowed the construction of consistent Minimal Left Right 
Supersymmetric models (MSLRMs) ITKt ITH] which were shown to naturally preserve R- 
parity and thus predict a stable LSP. Moreover these analyses brought clearly to the fore 
a theme that had been noticed from the beginningjH |T7j of the study of multiscale Susy 
GUTs : that generically in LR Susy GUTs there are light multiplets whose masses are 
protected by supersymmetry even though they are submultiplets of a GUT multiplet that 
breaks gauge symmetry at a high scale. As such they violate the conventional wisdom of 
the "survival principle" which estimates the masses of such submultiplets to be the same 
as the large vev. Thus it became clear [THj that Susy GUTs required the use of calculated 
not estimated masses for RG analyses. 

This ambling pace of theoretical development was forced by the epochal discovery of 
neutrino oscillations by Super-Kamiokande in 1997 and the rapid refinement of our knowl- 
edge of the parameters thereof ^H]- As is well known, the seesaw mass scale indicated by 
this discovery was in the range indicated by Grand Unification. Since LR models and GUTs 
containing them naturally incorporate the Seesaw mechanism, this gave a strong motiva- 
tion for taking up the detailed study of SO(IO) Susy Guts anew, particularly with the 
understanding provided by the natural class of MSLRMs developed earlier[T31 UHl HE]- A 
completely viable R-parity preserving Susy GUT based on the 10 © 126 FM (for Fermion 
Mass )Higgs system and an additional 45 © 54 © 126 AM (for Adjoint type Multiplet and 
euphony) system was then developed pOlj. The RG analysis carried out in this work already 
indicated that the use of calculated spectra would force together the various possible inter- 
mediate scales into a narrow range close to the GUT scale resulting in an effective "SU(5) 
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conspiracy". Moreover the problem of seeing how the various MSLRMs considered by us 
would fit into Susy GUTs had motivated a review of the various possibilities j^Il 1221 
this had again brought up the model based on the 210-plet as an important and interesting 
possibility. 

It thus became clear that a detailed calculation of the full GUT spectrum and couplings 
in various SO(IO) GUT models would be necessary if further progress was to be made and 
that this would require complete calculational control on the group theory of SO(IO) at the 
level needed by practical field theory calculations which was still lacking in spite of signal 
contributions [2ni 1211 • These techniques -based on an explicit decomposition of SO(IO) 
tensor and spinor labels into those of the maximal ("Pati-Salam") sub-group - were duly if 
laboriously developed by us|2Sl. At this time the model based on the 126 © 126 © 210 was 
again brought to the fore by us (211 and its old 13 I2Z1 claim to be called the minimal Susy 
GUT(MSGUT) was buttressed by an analysis of its parameter counting and the simplicity 
of its structure. 

The techniques developed by us [25] permitted first a partial calculation of the mass 
matrices of the MSSM doublets (1,2,±1) and triplets (3, 1,±|) [231 and then a com- 
plete calculation [2H1 of all the couplings and mass matrices of the MSGUT. With the 
same motivations two calculations, one in parallel and cross checking with ours j^j and 
another (201 quite separate, had commenced after |231. These calculations both used a some- 
what abstract inil method that permitted the calculation of the "clebsches" that entered 
the spectra of MSSM submultiplets of SO(IO) tensor (but not - so far - spinor) chiral 
supermultiplets (but not their couplings). After some controversy concerning the consis- 
tency of these three calculations |2S1 1211 IH21 (HHl 011 HEi a consensus seems to have emerged 
on their compatibility [^11 ESI '^'^^ withstanding notable differences in normalizations and 
phase conventions. In [23 I2H1 we also provided the complete (chiral and gauge) spectra, 
neutrino mass matrices, all gauge and chiral couplings and the effective d = 5 operators for 
Baryon violation in terms of GUT parameters. This laid the stage for a completely explicit 
RG based analysis of the MSGUT : an important example of which was the calculation of 
threshold effects based on these spectra[2Hl described in Section 3. 

Meanwhile the old theme of utilizing the very restricted structure of fermion Higgs 
couplings in the class of SO(IO) GUTs that used only 10 © 126 representations and renor- 
malizable couplings to make "predictions" concerning the neutrino mass sector was taken 
up again. It led [HI to a remarkably simple (2 generation) insight into the operation of 
the Type II mechanism that naturally generated large atmospheric neutrino mixing angles 
based on the observed approximate b — r unification in the MSSM extrapolated to Mx- 
This simple insight then provoked a detailed analysis of the fitting problem for dominant 
Type II seesaw mechanism and 3 generations[7 which was quite successful. A certain 
tension in the details of these fits was ameliorated by later workfTUl and disappeared |1U1 
when a 120-plet was also allowed to perturb the fit of the 10 + 126 slightly. Finally, very 
recently, in fact just in time for PLANCK05 Babu announced that they ^T] had succeeded 
in obtaining viable Type I, Type II and mixed fits to the (large mixing angle) neutrino 
mass data: which possibility had been neglected since the suggestion ofJH]. These sudden 
reversals are due to the extreme delicacy and complexity of the task of fitting accurately 
the masses of particles differing by a factor of upto 10^ in a multidimensional complex 
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parameter space where algorithms "shoot" obhviously past solution points. 

Obviously these successful fits of all low energy fermion "quadratic" data (extrapolated 
up to Mx ) in generic MSGUTs cry out to be married to a particular simple MSGUT 
in which all coefficients are specified and testable for viability. This is particularly so 
because only in a particular MSGUT, where the needed coefficients are fully specified PH] 
can one test the viability of a given FM fit (such as the Type II dominant onesjHl El HHj). 
Furthermore, as we shall see in detail below, the true output of the FM fit is really the 
embedding matrices that specify the relation between MSSM and MSGUT supermultiplets 
at the GUT scale. Given the phenomenologically expected quasi-identity of flavour mixing 
in the scalar and fermionic sector down to Electroweak scales, these matrices are the main 
missing ingredient for performing an informed and realistic calculation of Baryon decay and 
Lepton Flavour violation in the MSGUT. The combined cuts provided by gauge unification 
and fermion mass fit viability constrain the MSGUT and render it falsifiable . If it passed 
these tests with some region of its parameter space unscathed we would be in a position to 
actually specify the theory at the high scale: a task long thought to be impossible by those 
despondent at the continuing failure to detect proton decay. The key to this remarkable 
development is of course the remarkable uroborotic ( ovpo(3opo<, : the world-snake that 
swallows its tail) feature of the Seesaw which connects the physics of the ultralight and 
ultramassive particles. 

In Section II we review the structure and spontaneous symmetry breaking of the 
MSGUT and describe its decomposition to Pati-Salam labels. This is accompanied by an 
Appendix containing the complete spectrum of MSSM multiplets in the MSGUT -including 
gauge submultiplets. In Section III we describe the RG analysis of threshold effects and 
the type of cuts these corrections impose on the parameter space. In Section IV we review 
the fitting of Fermion masses and mixings with a view to supporting our assertion that the 
true output of this procedure is the specification of the embedding matrices. In Section V 
we plot the behaviour of the parameters that control Type I vs Type II dominance in the 
MSGUT and show that Type II dominance may be quite unlikely. We also use particular 
examples of solutions provided by one of the groups [10] that have performed the successful 
fitsjni Cni m] to plot the Type I and Type II contributions and again [HT] find indications 
that Type II is not easily dominant and that the overall size of neutrino masses tends to 
be too low except possibly in very special regions of parameter space. In the next section 
we transform our previously derived formula for the effective d = b Lagrangian for Baryon 
decay in the MSGUT to the MSSM using the embedding matrices. This operator is then 
ready for use to calculate Baryon decay. We also briefly comment on the relevance of our 
considerations to c? = 6 and Lepton Flavour violation. We conclude with a discussion of 
the outlook and the directions that will be pursued in the forthcoming detailed paper [¥1] 
on the application of the techniques indicated here to the pinning down of the MSGUT. 

2 Essentials of MSGUTs 

There are at present two main types of SO(IO) Susy GUT models namely renormalizable 
GUTs (RGUTs) [11121 120] of the classic type which invoke only gauge symmetry and preserve 
R-parity by maintaining the sharp distinction between Matter and Higgs Supermultiplets, 
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and non renormalizable GUTs (nRGUTs) which allow non renormalizable operators - 
particularly for generating fermion masses - but impose additional symmetries to perform 
the functions of R-parity and suppress unwanted behaviours allowed by the license granted 
by non renormalizability. Practitioners of the art are often sharply fixed in their preference 
for one type or the other and our own is obviously with the former. Many (talks by 
Babu,Malinsky,Mohapatra..) of the contributions on GUTs to PLANCK05 are of this type 
but the other type is also well represented (Pati,Raby...). 

Other distinguishing features of these two types of models are that RGUTs allow large 
representations such as the 126 Higgs to generate charged and neutral (via type I and 
Type II Seesaw) fermion masses but claim minimality on grounds of maximal economy 
of parameters. In contrast nRGUTs, allow only small Higgs representations and use non 
renormalizable operators (16|l6fj etc) to generate effective vevs in the 126, 120 channels. 
They achieve simplifications by invoking additional symmetries, whose rationale is however 
not always appreciated by non adherents. A very large number of works on nRGUTs [SQj 
have appeared and some of these have far advanced fitting programs : with claims to 
generate sample spectra with specific numbers for most (i.e ~ 10^) low energy MSSM 
parameters. Moreover they also seek to satisfy additional motivations like workable models 
of gauge mediated super symmetry breaking (HH]- 

The RGUTs, on the other hand are mostly at the stage of seeking a semi quantitative 
non- disagreement with the data as available at present. For example, since no superpartner 
has as yet actually been observed, and speculations on their possible masses range all the 
way from the classic 1 TeV up to almost the GUT scale, the value of the low energy param- 
eters renormalized up to the Susy breaking scale where the Susy RG equations take hold, 
can, realistically, be considered as known to at best an accuracy of 5 — 10%. Invoking the 
remarkable accuracy of experimental data at the scale Mz is of little avail since nothing is 
known of the size of the threshold corrections due to superpartners. It is this approach that 
we shall adopt in this paper and our quantitative analyses regarding quantities unprotected 
by any symmetry shall not lay claim to accuracy that may turn out to be quite spurious. 

A very important but controversial distinction between the two classes of theories is 
due to the fact that the one loop beta functions of the 126, 120 chiral supermultiplets 
are so large (35 and 28 respectively) that they alone overcome the contribution (-24) of 
the SO(IO) gauge supermultiplet and cause the gauge coupling to explode just above Mx 
the perturbative unification scale. By analyzing the complete two loop RG equations for 
the MSGUT above Mx we have shown P 5) that the large positive coefficient in the one 
loop gauge beta function precludes evasion of this difficulty by taking refuge in a weakly 
coupled fixed point before the gauge coupling explodes. Thus such theories effectively 
determine their own upper cutoff Kx ~ lO^^Gey . This feature is so inescapable that 
we proposed [Sni HOj to face it by making it a signal of a deeper non-perturbative feature 
of such Susy GUTs: namely that they dynamically break the GUT symmetry down to 
a smaller symmetry (such a.s H = SU{5) x U{1) or even H = G'123 ) at the scale Ax 
via a Supersymmetric condensation of gauginos in the coset SO{10)/H which drives a 
G123 preserving Chiral scalar condensate |40]. A new fundamental length A^^ below which 
SM probes cannot delve thus arises, endowing the particles of the SM with "hearts" i.e 
impenetrable cores with sizes ~ 10~^°cm. We have implemented [40 j the strongly coupled 
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Supersymmetric dynamics required by this type of scenario in a toy model based on SU{2) 
which can be easily generalized to at least the simple breaking 5*0(10) — > SU{b)xU{l)^^. 
Various fascinating perspectives including dual unification and induced gravity then beckon. 
However in view of the controversial nature of this proposal we shall not treat of it further 
here. This talk is aimed at indicating the feasibility of a promotion of the MSGUT to a 
falsifiable theory using the deep linkages between low and high scale physics inherent in 
the SO(IO) seesaw mechanism in combination with the traditional -relatively permissive- 
constraints of gauge unification. 

The MSGUT is the renormalizable globally supersymmetric 5*0(10) GUT whose chiral 
supermultiplets consist of "AM type" totally antisymmetric tensors : 210(<l'jjfci) 126(Sjjfc;m), 
12Q{Ylijkim){h i = 1...10) which break the GUT symmetry to the MSSM, together with 
Fermion mass (FM) Higgs lO-plet(Hj). The 126 plays a dual or AM-FM role since it also 
enables the generation of realistic charged fermion and neutrino masses and mixings (via 
the Type I and/or Type II mechanisms); three 16-plets ^'a(^ = 1, 2, 3) contain the matter 
including the three conjugate neutrinos (z^^). 

If in addition to the 10, 126 FM Higgs we also include a 120 -plet Higgs allowed 
by the SO(IO) multiplication rule 16 x 16 = 10 + 120 + 126, the GUT scale symmetry 
breaking is unchanged since the 120 contains no SM singlets. However it does contribute 
two additional pairs of (1, 2, ±1) doublets taking the MSSM Higgs doublet mass matrix TC 
from 4 X 4 to 6 X 6. The resulting theory is thus justly called the next to minimal Susy 
GUT (nMSGUT). 

The superpotential is 

W = -MhH^ + —^ijkl^ijkl + -^^^ijkl^klmn^mnij + -^^ijklm^ijklm 

TJ 1 

~(" ~^^^ijkl^ijmno^klmno ~^H{^jklm(^'y^ijklm ~l~ 'y'^ijklm'} 

The parameter counting is as follows : the 7 complex parameters m, M, Mh-, A, 77, 7 and 7 
can be relieved of 4 phases using the arbitrariness in the phases of the 4 SO(IO) multiplets 
210, 126, 126, 10 leaving 10 real parameters. The 12 complex Yukawa contained in the 
symmetric 3x3 matrices h\ f can be reduced to 15 real parameters by diagonalizing 
one linear combination of the two to a real diagonal form. In addition there is the gauge 
coupling. In all the MSGUT thus has exactly 26 non-soft parameters j^E]- Incidentally the 
MSSM also has 26 non-soft couplings. In[26j we have shown that this number of "hard" 
parameters is considerably less than any GUT attempting to perform the same tasks vis a 
vis fermion masses that the MSGUT accomplishes. The 'levity of the cognoscenti' provoked 
by the number 26 will recur again below ! Note however that one of these parameters, 
e.g Mhi is traded for the electroweak vev via the fine tuning condition that yields two 
light doublets and another Susy electro-weak parameter i.e tan/? must be taken as an 
additional given of the analysis since it is likely determined by dynamics dependent on the 
supersymmetry breaking parameters. 

The GUT scale vevs that break the gauge symmetry down to the SM symmetry (in the 
notation a, 6 = 1...6; a = 7.. .10 of |2S!) are P E] ((15, 1, l))2io : {4>abcd) = ^eabcdefeef, 
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((15, 1, 3))210 : (0aba/3) = ^^abe^p_ ] ((1, 1, 1))210 : {(PaPjs) = P^a^yS ! 

((10, 1, 3))y26 : (Si358o) = ^ ; ((10' 1' 3))i26 : (^24679) = ^- The vanishing of the D-terms 
of the SO(IO) gauge sector potential imposes only the condition \a\ = \a\. Except for the 
simpler cases corresponding to enhanced unbroken symmetry {SU{5) x U{1), SU{5) 
G3,2,2,B-L,G3,2,R,B-L ctc) j2Sl 121] this system of equations is essentially cubic and can be 
reduced to the single equationj^Hl for a variable x = — Acj/m, in terms of which the vevs 
a, uj, p, a, a are specified : 

8x^ - I5x^ + 14x - 3 = -^(1 - x f (2) 

where ^ = This exhibits the crucial importance of the parameter ^. 

Using the above vevs and the methods of[25| we calculated the complete gauge and 
chiral multiplet GUT scale spectra and couplings for the 52 different MSSM multiplet sets 
falling into 26 different MSSM multiplet types of which 18 are unmixed while the other 
8 types occur in multiple copies. {On the lighter note : the occurrences yet again of the 
'mystic' String Theory number 26 demonstrates that one can do just as well without string 
theory when searching for the fundamental theory !). These spectra may be found in the 
Appendix. 

Among the mass matrices exhibited is the all important 4x4 Higgs doublet mass 
matrix|23l2H] 'H. To keep one pair of these doublets light one tunes so that DetH = 0. 
This matrix can then be diagonalized by a bi-unitary transformation yielding thereby the 
coefficients describing the proportion of the doublet fields in 10, 126, 126, 210 present in 
the light doublets : which proportions are important for many phenomena. 



3 RG Analysis 

If the serendipity 132 113 113! of the MSSM gauge coupling unification at is to survive 
closer examination the MSGUT must answer the query : 

Are the one loop values of Siri^Ow o,nd Mx generically stable against superheavy thresh- 
old corrections ?. 

We follow the approach of Halljin] in which the mass of the lightest baryon number vi- 
olating superheavy gauge bosons is chosen as the common "physical" superheavy matching 
point {Mi = Mx) in the equations relating the MSSM couplings to the SO(IO) coupling |^ 

-IJF^ = L A + 8^b,ln^ + AnJ^^lnX, - 47rA,(Mx) (3) 
ai{Ms) aG{Mx) Ms ^ bj 

Seej'28j for details. In this approach, rather than enforcing unification at a point, it is 
recognized that above the scale Mx the effective theory changes to a Susy 50(10) model 
structured by the complex superheavy spectra which we have computed and which appears 
as unbroken SO(IO) only at the crudest resolution - here surpassed. Thus we compute 
the corrections to the three parameters LogioMx, sin'^9iY{Ms),aQ^{Mx) as a function of 
the MSGUT parameters and the answer to the question of stability of the perturbative 
unification is determined by the ranges of GUT parameters where these corrections are 
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consistent with the known or surmised data on LogioMx, sin^9w{Ms) and the consistency 
requirement that the SO(IO) theory remain perturbative after correction. We find the 
corrections 

M' 

A^^'^LogioMx) = .0217 + .0167 i5b[ + 3b'^ ~ 8b'^)Logioj^ (4) 

M' 

A.'^^^\sin^ew{,Ms)) = .00004 - .00024{4b[ - 9. Qb'^ + 5. Gb'^) Log w—^ (5) 

M' 

A(*^)(a^^(Mx)) = .1565 + .01832(56; + 362 + 1263)Lo^io—^ (6) 

Where b[ = IGn'^b'^ are 1-loop /? function coefficients {Pi = bigf) for multiplets with mass 
M' (a sum over representations is imphcit). 

These corrections are to be added to the one loop values corresponding to the successful 
gauge unification of the MSSM : 

Using the values 

a^aiMxy^ = 25.6 
a^\Ms) = 57.45 

the two loop corrections are 



= 10^^-^^ GeV ; Ms = 1 TeV 
a^^Ms) = 30.8 ; a^\Ms) = 11.04 (7) 



-.08 ; A'^-^"°P{sin^ewiMs)) = .0026 
-.546 (8) 

We see that in comparison with the large threshold effects to be expected in view of 
the number of heavy fields and their beta functions |46| 1^ the 2 loop corrections are quite 
small. 

A few remarks on the role of the parameters are in order. The parameter ^ = XM/rjm 
is the only numerical parameter that enters into the cubic eqn.Q that determines the 
parameter x in terms of which all the superheavy vevs are given. It is thus the most 
crucial determinant of the mass spectrum . The dependence of the threshold corrections on 
the parameters A, r], 7, 7 is comparatively mild except when coherent e.g when many masses 
are lowered together leading to ac explosion. The parameter ratio m/A can be extracted 
as the overall scale of the vevs. Since the threshold corrections we calculate are dependent 
only on (logarithms of) ratios of masses, the parameter m does not play any crucial role in 
our scan of the parameter space : it is simply fixed in terms of the mass My = Mx of the 
lightest superheavy vector particles mediating proton decay. 

With these formulae in hand we can explore the dependence of the threshold correc- 
tions on the "fast" parameter ^ in response to which the vevs and thus all the threshold 
corrections can gyrate wildly (see plots for the real solution for real ^ below ), and the 
"slow diagonal " parameters A, 77 whose lowering tends to make fields light and thus give 
large negative corrections to a^^. In addition there are the "slow off diagonal" parameters 
7,7 whose effect seems quite mild. 



M 



X 



^ Ms 
A^-^""Paa\Mx] 
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Sin^Ow is now very accurately known jIH] at Mz'- s| = .23120 ± .00015 or an error of 
less than .065%. Similarly the value = 127.906 ± .0019 has an error of only .015%. 
The main uncertainty in the data (besides Mh ) at Mz is in a^^Mz) : ~ 1.5%. However 
the same is obviously not true of the overall sparticle mass scale or the intra sparticle mass 
splittings : with values in current speculation ranging from ~ 1 TeV to 10^° TeV ! Thus at 
least until the first superpartners are observed the assumption that we know the values of 
the ai{Ms) with anything like the precision at Mz is quite unjustified. A rough estimate 
of the uncertainty in Sin^OwiMs) gives numbers in the ball park of 1% to 10% and we 
shall not pretend to have any license to impose stronger constraints on our parameters. For 
the uncertainty in LoqiqMx however there is a quite stringent Ms independent constraint 
coming from the requirement that gauge mediated proton decay in channels like p —>■ c^tt^ 
obey the current bounds Tp^g+TrO > lO^^yr. Thus ALoqiqMx > —1 is a very firm constraint 
that we may rely upon. Finally it is clear that the applicability of the perturbative formulae 
used in our treatment will need detailed scrutiny if the fractional change in aG{Mx)~^ is 
greater than about 25%(particularly if the change is neg ative). So \ AaG{Mx)~^\ < lO(say) 
is a limitation we tentatively observe as a limit on the range of validity of our calculation 
( rather than the theory itself : since a detailed examination may permit one to handle 
larger changes - particularly positive ones - in aciMx)"^)- In view of these relatively loose 
bounds on the acceptable changes associated with threshold effects we cannot expect that 
gauge constraints can limit the allowed parameter space very precisely. Nevertheless, these 
constraints used in conjunction with additional very significant information arising from 
the fitting of Fermion masses and (in more model dependent ways) d = 5 Baryon violation 
and Lepton flavour violation could allow us to obtain a complete semi-quantitative picture 
of viable regions, if any, of the MSGUT parameter space. 

We now present examples showing how such a semi-quantitative mapping of the RG 
constrained topography of the GUT scale parameter space may be obtained using the 
formulae given above. An attempt at an exhaustive mapping would be premature till the 
recent available FM fits fOl have been digested. Later we shall present our program 
of further constraining the parameter space using the constraints of the fit[Hl IHl CHI lUI of 
Fermion masses and also due to the consequent specification of Baryon decay operators. 
It should be kept in mind that at any ^ the cubic equation (j21 ) has three solutions any of 
which is in principle exploitable for defining a vacuum. 

Consider first the plots of the threshold changes A'^'^'^\aQ^{Mx)) and A^*'^\LogioMx) 
vs ^ after (arbitrarily) setting {A, r], 7, 7} = {.7, .5, .3, .2} which are shown as Fig. 1 for the 
real solution of the eqn© for real ^. 

The "twin towers" due to singularities at ,^ = —5, 10 in the plot of adMx)'^ arise 
from the SU{5) invariant real vevs at these values of ^. Similarly the negative spike at 
^ = —2/3 corresponds to a real solution with SU{5) jupped x f/(l) symmetryj^Hl I^H]- As 
emphasized by us in[28j the plots of the threshold corrections show very sharply defined 
features corresponding to the special behaviour near points of enhanced symmetry and 
thus provide a rapid way of scanning the topography of the parameter space. For these 
moderately large values of A,// the lower cut at AaG{Mx)~^ = —10 essentially allows all 
possible values of ^ though the upper cut does advise caution around the special symmetry 
points with SU{5) symmetry. However the graph of ALogioMx with a cut at ALogioMx = 
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— 1 immediately rules out most of the region ^ G (.25, 8.6) due to the ultrarapid gauge boson 
mediated Proton decay in that range of ^. 

When we lower the value of A or many particles became light so that the entire 
aG{Mx)~^ vs ^ plot shifts downward. This is seen in Fig. 2 where we repeat the same 
Plot as Fig 1. but now with A = .1. As a result larger C, values are allowed on grounds 
of limiting the change in aG{Mx)~^ ■ On the other hand the behaviour of ALogioMx in 
response to a decrease in the diagonal slow parameters A, 77 is quite mild. Thus although 
lowering A does tend to make ALogioMx less negative in ^ > region, the large positive 
change in adMx) would rule out the small A, 77 region. The condition ALogioMx > — 1 
required by proton stability rules out the region .25 < ^ < 8.6 for the real case (see Figs. 
l.,2.). 

Next we plot ASin'^6w{Ms) versus for moderate and very small A (Fig 3.). The 
change in A has practically no effect. We also see that there are ranges of C, where the 
change in Sin'^9w{Ms) is less that 10% and that these ranges are only slightly affected by 
the change in A. However the region 2.8 < ^ < 8 which was excluded by ALogioMx > — 1 
is also excluded by the large change ASin'^9w{Ms). 

The further evolution of A(*^)(a(.^(Mx)) and A(*^)(Lo^ioMx) as A is decreased to .01 
is similar as can be seen in Fig. 4. Lowering rj has effects very similar to those of lowering 
A since the effect of smaller values for these diagonal couplings is to lower the masses of 
whole sets of multiplets and therefore raise adMx) sharply. The effects of the off diagonal 
couplings 7 and 7 are much milder. 

With ^ real there are also two complex (mutually conjugate) solutions of the basic cubic 
equation 0. Examples of these are displayed as Fig. 5,6 for moderate and small A. We 
observe that the behaviour of the complex solution is smoother than the real one and that 
apart from the spikes observed at = —5 there is hardly any sharp feature to be seen. The 
effect of decreasing A or 77 on A^^^\aQ^ (Mx)) is however even stronger than in the case of 
real x as may be seen in Fig 6, thus restricting the viable magnitudes of A, rj to moderately 
large values again. The corrections to Sin^9^/{Ms) for the complex solution are shown 
for moderate and very small values of A as Fig. 7. One can also consider complex values 
of the parameter ^ and the three solutions in that case. The behaviour is quite similar 
to that shown for the complex solutions for real ^ but needs detailed and comprehensive 
examinat ion |lTj . 

We have exhibited these graphs to give a sense of the structure visible once one ramps up 
the resolution of analysis to reveal the finestructure hidden within the bland impressiveness 
of "supersymmetric unification at a point". As already noted long since |37j the ambiguities 
associated with superheavy thresholds would not allow one to predict the effective scale 
of superpartner masses or the unification scale even if the low energy values ai{Mz) were 
known exactly. In fact, following Hall[45j we have chosen not to treat the unification be- 
yond leading order by imposing unification at a point but rather in terms of quantifying 
the ambiguities in a^^(Mx), Sin^OwiMs) and LogioMx caused by the finestructure of 
unification scale mass spectra. The range of behaviours exhibited make it unlikely that 
the constraints of gauge unification alone will rule out this minimal Susy SO(IO) GUT or 
fix its parameters. However when taken together with the rest of the low energy data, the 
MS GUT provides a well defined and calculable framework within which significant ques- 
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tions regarding the viable GUT scale structures can be posed and answered. Furthermore 
the processjHl El CDl HI] of fitting the highly structured fermion data to the relatively few 
parameters of the MSGUT excavates crucially important information regarding the em- 
bedding of the MSSM within the MSGUT. This information (concerning mixing angles 
of various sorts) is critical for determining the precise predictions of the MSGUT for both 
d = 5,6 baryon number violating operators in the effective MSSM as well as the predictions 
for Lepton Flavour violation. 



4 FM Fitting Frenzy 

We have already reviewed the sequence of developments preceding the current focus of 
interest on the fitting of fermion mass and mixing data in the MSGUTjHl 13 |H1 UHl ITT| . 
The fitting program itself has used only the form of the fermion mass formulae in the 
MSGUT (which follows from the use of only the 10 + 126 representations) rather than 
the specific formulae for the coefficients in the fermion masses dictated by the MSGUT 
superpotential. Our concern here is not with the actual values of the successful fits but 
rather their implications when combined with the structure of the MSGUT. We therefore 
review the fitting procedure from our particular viewpoint. To begin with the "Clebsch 
" coefficients for the couplings of the 16 x 16 SO(IO) chiral spinors to the 10, 126 irreps 
were calculated as a part of our explicit decomposition of 5*0(10) in terms of Pati-Salam 
labels : One obtains : 

= -2\/2h'^j^hi [(IaQb + gaLb] + 2 V2/i'^^/ii [uaQb + i^aLb] 

+ (9) 



= 4a/2/^^[-^{/12(c/a<5b - 'itA^B) - h2{uAQB " SuaLb)} 

+ A/'^bI-'^iG^uai^b) + V20LaLb) + 



(10) 

where the alphabetical naming convention regarding the subcomponents of the Higgs 
multiplets is given in the Appendix and in detail in[2H]- From the properties of the SO(IO) 
Clifford algebra it follows that the Yukawa coupling matrices H'ab, /'ab^^j B = 1,2, 3) are 
symmetric {h' ab = h' ba-, J'ab — f'sA) complex matrices. Therefore the freedom to make 
unitary changes of basis allows one to chose one linear combination of the matrices , f to 
be diagonal. We shall choose the basis where /' is diagonal since it proves convenient when 
analyzing the Seesaw mechanism but our conclusions are independent of any such choice. 
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To obtain the formulae for the charged fermion masses from the above decomposition j2H] 
one first needs to define the (G32i(l, 2, ±1)) multiplets H^^\ H^^^ which are the (hght) 
MSSM Higgs doublet pair. This is achieved by imposing the condition 



Detn ~ O(Mvk) 

on the doublet mass matrix Ti which occurs in the quadratic terms of the superpotential 
when expanding around the superheavy vevs: W = h7ih+.... This amounts to a fine tuning 
of(say) the mass parameter Mh of the 10-plet Higgs. 

The 4x4 matrix can be diagonalized by a bi-unitary transformation : j^Hl 1211 

1231211 from the 4 pairs of Higgs doublets h^^\ h^'^^ arising from the SO(IO) fields to a new 
set H^^\ H^^^ of fields in terms of which the doublet mass terms are diagonal. 



U^nU = Diag{ra'^\7n§,--) 

/i« = UijH^^'^ ; h^^^ = UijH^^^ (11) 

Then it is clear that in the effective theory at low energies the GUT Higgs doublets 
h^^\h^^^ are present in in a proportion determined by the first columns of the 

matrices f/, U : 



E«Mx : h^'^ aiH^^^ ; ai = Un 

/i« ^ a,iy« ; a, = Ua (12) 

Thus the formulae for the charged fermion masses in the MSGUT are|2H] : 

M'^ = rji + r2f 

M' = ri/i-Sra/ (13) 
= h + f 



where 

h = 2V2h'vaismf] 

f = -A^tf'a2 sin P (14) 

ri = — cot p 

•^2 , n 
r2 = — cot p 

here {h', /') are the couplings in the MSGUT superpotential. 

Similarly the Majorana mass of the 'right handed neutrinos ' i.e of the fields u = u'^ is 
read off from the decomposition given above j28j : 

= -AtV2f'^^ < Si^^ >= AV2f^^a (15) 
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and is of the order of the Unification scale or somewhat less. The left handed or SM 
neutrinos receive a direct Majorana mass from the so called Type II seesaw mechanism|49j 
when the left handed triplet O contained in the 126 field obtains a vev. One obtains [22] : 

= ^V2f'^B < O'' >= Sif^B <0-> (16) 

The vev < 0_ > : (10,3^,1)2 arises from a tadpole following SU{2)l breaking (see 
below) . 

The final component of the Seesaw is the Neutrino Dirac mass which links the left and 
right handed neutrinos : 
Dirac mass : 

m^^ = 2V2h'j,j, < h^^^ > +4iVQf'AB < > (17) 



To determine the Majorana mass terms of the left handed neutrinos in the effective 
MSSM we must eliminate the superheavy u and evaluate the Type II Seesaw tadpole. One 
then obtains |2H! (some factors of \/2 have been corrected relative to eqns(77-79) of^Hj). 

Ml = -]m^Mr^m^ = -u{h - 3f)r\h - 3/) = -r,h (18) 
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Mi' = ^ifAB < O- >= r,f (19) 



ra = -2tV3{ai-f + 2V37]a2){—){^)Stnf3 

a2 Mo 

-ia2Sin(3 v 

= -= (20) 

4^3 a 

Mo = 2{M + r]{3a - p))) 

Thus we see that the fermion mass formulae are completely determined in terms of the GUT 
scale breaking parameters(^, A, ?7, 7, 7, m), the 10 + 126 Yukawa couplings(15 parameters) 
and the low energy parameters vew = 174GeV, tanp. 

To perform the fit, we must match the fermion masses and mixings of the MSSM RG- 
extrapolated to the GUT scale Mx '■ 

Lm%m{Mx) = r^Dil + u''^DuU + (r^Dd 

+ I ^v + U ^d' + u^{VD^V^)iy + ■■■ (21) 

with the effective theory derived from the MSGUT by integrating out the heavy fields at 
Mx- Here the Di^u,d,u-, C, V : are the diagonal fermion mass matrices (with mass eigenvalue 
components /^i, ma, c^a, ^a), the CKM mixing matrix and the PMNS matrix (Majorana 
Neutrino mixing) at the scale Mx in some fixed phase convention for the MSSM masses 
and mixings (either at Mz or at Mx) and d' = Cd. For example the convention that 
all the diagonal masses are real and positive. Currently each group of FM fitters uses 
idiosyncratic phase conventions. A standard format presentation of the data of the MSSM 
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at Mx is needed and is being pursued. The corresponding quadratic terms in the effective 
tlieory derived from tlie GUT are (GUT fields carry carets) : 



r(2) 
^GUT 



wliere 



= {h + f)u + (f{hri + fr2)d + l- (hn - 3r2f)l 



Mu = ra/ - r^{h + f)r\h + f) = nf - rS 



(22) 



(23) 



Wfien equating tlie two quadratic forms we must allow for Unitary transformations 
between the fields in the two Lagrangians : 



-(2) 
^GUT 



'MSSM 



{Mx) 




= Q 



( u 




u 



The unitary matrices Q, £, V(uc), V(dc), V(ic) describe the embedding of the extrapo- 
lated MSSM within the MSGUT. The {d = 5,6) effective lagrangian ( Cff/\il;)) must 
be transformed to the extrapolated MSSM basis in order to derive rates for such exotic 
processes. Thus it is clear that these matrices are neither unphysical nor conventional once 
the conventions of the MSSM parameters are fixed. In fact we argue that the crucial infor- 
mation given to us by the fitting procedure is not a prediction of the neutrino masses and 
mixings but rather information on these embedding matrices. 

Since the 10, 126 Yukawas are symmetric it follows that only two (say Q, £ ) of these 
matrices are independent while the others (V(„c), V(dc), V(ic)) can be determined in terms 
of the two chosen to be independent plus arbitrary diagonal unitary matrices 
Once this is done we obtain : 



hr,+r2f = VJcL'<iCtQ = 7^^L>'7^ = 7^^$,L>,7^ 



hri-3r2f = Y^DiC ^ D[C = C^^iDiC 



(24) 



r^f-rS = C^VD.V^C 
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where we have defined 1Z = C^Q. The phase freedoms $„, $d have been found P ITT] 
to be important for arranging the tunings that underhe the successful Type I and Type II 
fits of the fermion masses. However the phases play no important role so far since the 
phases in the PMNS matrix V are unknown at present. If we reabsorb the phases in the 
equation from Mi in L then we also need to redefine the PMNS matrix V to reabsorb them 
in the neutrino mass equation (last of equations (j^ V This ambiguity should be kept in 
mind when deducing predictions of Leptonic phases from the fit , but does not play any role 
at this stage. 

We remind the reader that all the parameters rj; i = 1,2,3,4 are known in terms of 
the MSGUT parameters but the explicit form is not used or invoked when solving the 
fitting problemjnilH!- Rather the parameters and the Yukawa couplings are determined 
in terms of the extrapolated experimental data. Thus the compatibility of the FM fits 
with the MSGUT remains to be verified for each prima facie viable fit since the ability 
of the MSGUT to reach the required values of the rj simultaneously while preserving the 
constraints of perturbative gauge unification is not obvious. 

4.1 Solution of the fitting problem 

The equations for the down fermion and charged lepton masses may be immediately solved 
to yield 

(25) 

we define 



V = n*C^, u = C^D'^C (26) 

for later convenience. 

Consider first the casejHl QHl where the mixing matrices, (hatted) Yukawa couplings 
and parameters are all assumed real. Moreover the phases $d are simply signs and the 
CKM phase 5 is also a sign which in fact is found to be minus in "successful" fits. 
Then solving the Up quark equations and eliminating ri,r2 using 23 component and trace 
we can put it in the form 

X^Xi^ VDiV^ = ^iu- ^D',) + ^D', (27) 

"23 J-d J-d 

Here T'^ = tr[D'j]. Notice that the non diagonality of Xi is completely driven by the 
non-diagonality of the matrix u' which in turn follows from that of the CKM matrix. 
FoUowingjU] it is convenient to rescale each of the diagonal fermion mass matrices by the 
mass of the third generation fermion of that type to get a dimensionless (tilde-ed) form of 
the equations : 
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Df = mf.Df- 



(28) 



X, ^ VDp^ = ^(^ - ^D^,) + ^D^, (29) 

■"23 -^d' -'J' 

The matrix X^- has eigenvalues /i,2 and /a = 1 by construction and hence it follows that 
we must have 

det(Xi - ~kh) = (30) 

Thus we obtain three coupled non-linear equations for ^"23 and two other quantities 
which are conveniently chosen to be di^2- We can solve numerically for X23 and di^2, given 
any set of up quark and charged lepton masses together with CKM data. Obviously only 
solutions within the error bars(usually allowed to be l-o") for the down quark masses are 
accepted. Note however that inasmuch as there is a strong inter-generational hierarchy for 
the charged fermion masses the numerical solution of the three coupled non-linear equations 
is an extremely delicate operation requiring utmost care and diligence to find the correct 
solutions. Since brute shooting usually fails in multidimensional problems of this delicacy 
approximate analytic solutions obtained by expanding in the light fermion masses are used 
to guide the numerical search for solutions. 

Assuming this has been done the matrix Xi is completely determined so that on diag- 
onalizing it one obtains the crucial matrix T> — TZC^. With V in hand one can proceed by 
choosing the convenient /-diagonal basis mentioned earlier. Rescaling 

j^f/mt^f^C^X^C (31) 
where we have defined Xj and unitary !F by 

Xf = V^D'p - Dj{^) = TDjT'^ (32) 

Since / is diagonal by choice of basis it immediately follows that 

T ^n^VC (33) 

Now since C, TZ are known and since fi_2 = ^^1,2 are calculable in terms of Xf it is 
clear that we have also determined the Yukawa coupling h : 



h = h/mt = riC^[W^D'p + Di{—)]C (34) 

In other words one finds that the fitting of charged fermion mass ratios requires tuning 
of the down quark mass ratios di^2 to less than one part in 10^^ for given precise values of 
up quark and charged lepton masses together with CKM data and yields the dimensionless 

matrices /i, / and the exotic embedding matrices TZ and, given mj, /i, /. 

In the realistic case when the parameters are complex a similar but numerically even 
more difficult procedure is followed. Like the sign freedom in the real case the phase freedom 
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of choosing $rf is found P CDl HI] to be crucial to obtaining a successful fit. The u 
equation still has the same form 

Xi = VDiV'^ = pu + qDd (35) 

but now TrXi ^ TrDi. 

To proceed we solve the 23, 33 components to eliminate p, q and get 

Xi = VDp^ = ^^ + (X33-i^^33)i5j (36) 

U23 U23 

Then to determine V we must diagonalize 

XjXj = VDfV^ (37) 

This matrix has eigenvalues 1 and an inspection of its explicit form shows that 
it requires knowledge of IX23I, |X33|,0x = Arg{X23) — Arg^Xs:^) and a choice PI ITUl ITT] of 
the quark mass phases The numerical results obtained by these authors are then 

a specification of the right hand side of eqn (jHH|) consistent with some acceptable values of 
the charged lepton masses. We refer the reader to the original papers for the procedure 
for fixing the unknown phases. Once this rather horrendous numerical problem has been 
solved (its trickiness accounts for the fact that even 12 years after it's solution was first 
explicitly attemptedlHj this system continues to throw out surprises jTT]) one can proceed 
essentially as before : diagonalizing the rhs of eqn ()36|) numerically thus yields the matrix 
V. 

Since / is diagonal and real by convention one writes 

f=~f/mt = f2C^XjC (38) 

where 

777 

Xj = V^D^,V* - Dii-^) = J^DjT^ (39) 
Since Xj is symmetric , JF is Unitary and -Dj is real like /, it follows that 

I=\h\Dj (40) 

so that 

C = re — ^ (41) 

where T is found by diagonalizing ^j^j = J^DjJ^^ 
Finally 

n = V*C ; Q = Cn (42) 

In summary : for given up and charged lepton masses and given CKM mixing angles, by 
tuning di^2i ^ckm within their allowed ranges and for a certain choice of the phases $„, $d, 
one completely fits the charged fermion masses and determines the FM Yukawa couplings 
h' , f of the GUT along with the exotic embedding matrices £, Q, 
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4.2 Fitting Neutrino Masses 

The current FM Fitting Furore jHl El was triggered by the remarkably simple obser- 
vation of[H] regarding naturalness of large atmospheric mixing angles in MSGUTs with 
dominant Type II seesaw mechanisms given the near equality of the r lepton and bottom 
quark masses at Mx- If one assumes Type II domination i.e rs >> r4 the Seesaw formula 
simplifies to just 

f = (^—)C^VD^V^C (43) 
Since we chose a/-diagonal basis it immediately follows that 

V = Ce^^ ; D, = \rs\f (44) 

Thus the neutrino mixing angles and ratio of mass squared splittings can be determined 
under these assumptions since we know f,C. 

In the general case the neutrino mass fitting equations take the scaled form : 

ATDj^AT^ = /- C-^)n = {J—)C^VD,V^C (45) 

where 

n = n/mt = {1 - ?>I)T\^ " 3?) (46) 
Thus the mixing matrix and neutrino masses are also completely determined : 

iArgirri ) 

V = C*Me^^ ■ D, = mt\r^\D^ (47) 

These fitting problems have been formulated and solved with increasing refinement by 
a number of authors|ni lEl E]- It seemsfTI] that Type II dominant as well as Type I and 
Type I plus Type II combined solutions (which are not perturbations of Type II dominant 
solutions ) can be found in the complex case. Further solutions probably still remain to 
be found and the possible solutions definitely still need to be compactly characterized and 
parameterized. The fits achieved so far already motivate a detailed examination of what 
type of Seesaw is actually allowed by the MSGUT in various regions of its parameter space. 
We have argued that in each case the exotic embedding matrices are -like the 'philosophers 
stone' - the so far unregarded true product of the fitting calculation. We emphasize that 
in practice one accepts solutions which give v oscillation parameters which lie within the 
1 — 0" ranges around the central values of the fermion mass and mixing parameters so that 
the true output of the fitting calculation are the previously completely unknown embedding 
matrices Q, £, $d which specify how an MSSM(Mx) with given conventions lies within 
the MSGUT. These matrices are crucial for pinning down the prediction of /S.B ^ pro- 
cesses in the MSGUT. Before considering that aspect however we turn to the question of 
what kinds of solutions are compatible with the MSGUT. 
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5 Scanning the MSGUT for Neutrino Masses and Mix- 
ings 

The possibility of large PMNS mixing angles is well understood in the Type II dominant 
case jHl EI, where it appears as a natural corollary of the approximate unification of the 
running bottom quark and tau lepton masses at scales 0{Mx) due to the 3 fold faster 
evolution of the bottom quark mass. As is evident from eqn. (j^3j) the parameter which 
controls the strength of Type I versus Type II Seesaw in the MSGUT is the ratio of the 
coefficients and r^. Complete Type II dominance requires << r^. To illustrate how 
the MSGUT yields information on this ratio we work with an example of a quasi-realistic 
real Type II fit that ignores the CP violating phase, which was kindly provided to us 
by S.Bertolini and M. Malinsky. Although Type II and only semi-realistic it should be 
emphasized that the values of h, f given will be rather typical since they are fixed by the 
charged fermion mass data and further selection is then imposed based on compatibility 
with the neutrino data. 

The data of the example solution ( at Mx and for tan /3 = 10 ) is 

Du = {0.785556, -191.546, 70000}Me\/ 

A = {0.3585,75.7434, 1290.8}Me\/ 

c/3 = mb = 1138.07Me1/ 

{SinO'l^, Sine^23, SinO^^^} = {0.2229,0.03652,-0.00319} (48) 

Notice the negative sign of U3 and 9^^ (this corresponds to taking 6ckm = tt), moreover 
the fitting procedure gives 

X23 = -0.14325 ; di = -.001105 ; = -.02747 (49) 
so that the sign ambiguity $d is also fixed. 

/ 0.98953 -0.136941 0.045582 \ 
V = -0.128508 -0.979735 -0.153642 (50) 
V -0.0656981 -0.146176 0.987075 / 

Supplying = 1138.07 MeV allows one to calculate 

/ 0.794035 -0.582472 0.173881 \ 
£ = = -0.563923 -0.599057 0.568437 (51) 
V 0.226934 0.549415 0.804142 / 

Then the matrices TZ, Q follow : 



/ 0.87329 
0.415588 

V 0.254266 

/ 0.943138 
0.219729 

^ 0.24942 



-0.469294 
0.577356 
0.668149 

-0.330924 
0.691351 
0.642279 



0.130873 \ 
-0.702813 

0.699233 / 
-0.0313091 \ 

-0.688297 

0.724753 / 
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The scaled Yukawa couplings are 



/ 



/ -0.000595629 ° \ 

0.0020326 

V -0.00804198 / 

/ 0.0626836 0.159778 0.181181 

0.159778 0.409183 0.466796 | (53) 
\ 0.181181 0.466796 0.532013 



The matrix n multiplying — r4 in the Type I mass is 



n 



( 1.4996 3.8747 4.55331 

3.8747 9.98038 11.6875 | (54) 
V 4.55331 11.6875 13.63 



From these the Type II mixing angles and the ratio of the 23 and 12 sector mass squared 
splittings is found to be 

{Sin^{2e(^),Sin^{2e^^),Sine(^} = {0.90981,0.88873,0.17387} 

Aml^/Amls = 0.06237 (55) 

The Yukawa couplings obtained are typical of Type II fits. It is clear that unless r^/r^ 
is very small the Type I term will dominate completely. To see just how small should 
be we plot the neutrino mixing angles SirP26i2, Sin^2923 versus this ratio as Fig 8. 

From Fig. 8 it is evident that a value of |r4/r3| > .0003 causes a collapse of the large 
mixing angles of the Type II dominant solution. Thus we should not expect a Type II 
solution to work in the MSGUT unless R = r^/r^ (which is completely determined by the 
GUT as given in the previous section ) obeys \R\ << 10^^. The MSGUT formulae for 
r3, r4 are given in eqn ()21|) 

Using these we plot the ratio R versus ^ and find that its typical value is ~ 10^^ or 
more, not 10~^ or less. Illustrative plots for the real and complex solutions of the cubic 
eqn.(j2I), for real ^ and A = .7 are shown as Fig 9. 

From the plots it is evident that in the real case R has a chance of being of the required 
small size only in the window near ^ = —.7 and between about ^ = 1 and ^ = 4. However 
in the former region one finds that the corrections to Siri^Oy/ diverge, while the region 
^ G (1,4) is not allowed by the requirement that ALogioMx > -1- On the other hand, 
in the complex case, the ratio seems bounded below by 1 ! Thus we see that if Type II 
fits were the only allowed ones then the combination of the gauge unification requirements 
with those imposed by neutrino mass phenomenology tend to rule out the MSGUT based 
on this class of solutions. This conclusion however still requires a more thorough study 
of all possibilities. Moreover recent work ^T] has shown that in fact large mixing angles 
can be achieved even in Type I solutions and Type I-II combined solutions (which are far 
from pure Type II solutions). Our intent, in this talk, is not to provide an exhaustive 
survey of the possibilities but only to illustrate how the combined requirements of neutrino 
oscillations and baryon stability can severely constrain the MSGUT over its parameter 
space. A detailed survey of the MSGUT for each of the three solutions of the basic cubic 
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equation to find in which regions, if any, the ratio R matches that required by the FM fit 
is quite feasible using our methods, and will be reported elsewhere ^Tj. 

Note that while the FM fits described in the previous section do not determine the over 
all mass scale of the neutrinos since the input data does not contain this information, the 

same is not true in the context of the MSGUT. Using the dimensionless versions / and h 
given by the real Type II FM fitting analysis as a guide to typical values of the Yukawa 
couplings allows one to compute the magnitudes of all three neutrino masses for each type 
of fit. When this is done another problem becomes apparent : except possibly for narrow 
ranges of ^ the largest neutrino mass (Type I or II) is much smaller than the mass splitting 
known from atmospheric neutrino oscillations. Thus even if one can find a region where 
the ratio of mass splittings and mixing angles for neutrinos are in the allowed region the 
additional consistency constraint that \{m^)max\ > \{^^)atmos\ alone can exclude most of 
the parameter space ! An illustrative plot of the maximum type I and Type II neutrino 
masses for the Yukawa coupling matrices that arises in the Bertolini and Malinsky solution 
used by us for illustration is given for the real solution as Fig. 10. Due to the cut at ^ = 8.6 
imposed by the condition ALoqiqMx > — 1 we see that there is no region with large enough 
values of \{mu)max\- Similarly in the complex case we get the plots shown in Fig. 11. As 
expected the Type I dominates completely. However the mass magnitudes for Type I tend 
to be too small : as can be seen from Fig. 11 and the magnifications shown in Fig 12. 
Since the shortfall is only a factor of 10 and we have not yet used the complex Type I fits 
found in^J it would be premature to rule out Type I fits. Nevertheless a certain tension 
is apparent. 

If the impossibility of large enough overall neutrino mass scale is borne out by a com- 
prehensive analysislHj it would require a revision of the perturbative MSGUT. Similar 
arguments were used in[SZj (for the Type II FM fit case only) to motivate an extension of 
the model by introducing an addition 54 plet so as raise the value of the 126 vev to a high 
scale to allow Type II to dominate and yield a large enough neutrino mass. Note however 
that in that case the minimality of the MSGUT is seriously diluted necessitating a complete 
reanalysis which has not been performed so far. An alternative to this extension may be 
to implement a non perturbative mechanism |40l 1^ based on dynamical GUT symmetry 
breaking down to SU{5) x f/(l) at Ax ~ lO^'^GeV^ and then use a Type I fit with small 
126 vev (T to achieve a larger value of the overall neutrino mass scale. 

6 AB 7^ IN MSGUT 

Finally we briefly indicate the d = 5 baryon decay operators determined by the FM fit. 
If any FM fit is found to be consistent with all the constraints discussed above then the 
computation of the actual Baryon Decay predictions will become a worthwhile exercise. 

The basis for this is the effective superpotential that arises when the superheavy Higgs 
triplets are integrated out from the theory ^K] ESj 

= ^ABCDi^f^QAQsQcLD) + RABCDieeAfuB^cdn) (56) 
Labcd = S-^hABhcD + S-^hABfcD + <S2^ fABhcD + Jab fen (57) 
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Rabcd — S^^hABhcD — S^hABfcD — fABhcD + fAsfcD 

- lV2S^^fABhcD + lV2S2''fABfcD (58) 

Here S = T^^ where T is the mass matrix in the [3, 1, ±2/3]-sector which is the repre- 
sentation type that mediates ci = 5 baryon decay. The Yukawa coefficients hAB,fAB are 
related to those in the superpotential by Hab = Sy^/i'^^, /ab = 4v^/^^ 

Substituting for the GUT superfields (with carets) in terms of the embedding matrices 
Q, 7^, C determined by the FM fit 



Ql = Q^Ql Ll = C^Ll (59) 

u = uc = Q^Kuc d = 4 = 7^t$rf4 i = L = c^'^*iic 

one obtains the coefficients of mass eigenstate MSSM fields in the AB ^ superpotential 
to be 

Labcd = QaeQ^bfQcg^^dh^efgh 

Rabcd = C\EQ^j,j,Qla{'R-^)DHmEE{K)FF{K)GGm)HHREFGH (60) 

A similar calculation[41^ can be done to determine the effective operators for d = 6 
baryon violating operators, which will be relevant even if the Supersymmetry breaking 
scale eventually turns out to be large enough (> lOOTeV) to exclude any observable d = 5 
baryon decay. Similar care needs to be exercised when studying lepton flavour violation. 



7 Conclusions and Outlook 

The MSGUT based on the 210 © 126 © 126 © 10 Higgs system is the simplest Supersym- 
metric GUT that elegantly realizes the classic program of Grand Unification and which is 
prima facie compatible with all known data. Its symmetry breaking structure is so simple 
as to permit an explicit analysis of its mass spectrum at the GUT scale and an evaluation 
therefrom of the threshold corrections and mixing matrices relevant to various phenomeno- 
logically important quantities. Since it has the least number of parameters of any theory 
that accomplishes as much this theory currently merits the name of the minimal supersym- 
metric GUT or MSGUT. The same simplicity and analyzability of GUT scale structure 
also applies to the theory with an additional 120-plet, since it contains no standard model 
singlets, and thus justifies calling it the nMSGUT. 

The small number of Yukawa couplings of the MSGUT makes the fit to the now well 
characterized fermion mass spectra very tight so that the combined constraints of this ffi 
and the preservation of the gauge unification observed at one loop in the MSSM may well be 
enough to rule out most or even all of the parameter space of this theory. A corresponding 
investigation is also possible for the nMSGUT and both are in progress |41j. 

We have emphasized the need for clarity regarding the flavour basis used when perform- 
ing the FM ffi. When this is maintained it becomes obvious that the true outputs of the 
FM fitting calculation are really the embedding matrices that define how the MSSM, with 
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fixed phase conventions, lies within the MSGUT. With these matrices determined by the 
FM fit in hand, one will be in a position to perform a much more reliable calculation of 
d = 5,6 Baryon decay in the MSGUT with the Susy breaking scale as the chief remaining 
uncertainty. This realization underlines and emphasizes the organic connection between 
the physics of neutrino mass and Baryon violationjSU] i.e between the physics of ultra- 
heavy and ultra-light particles which is the most intriguing implication of the discoveries 
of Super-Kamiokande. 

Our discussion concerning embedding angles also has implications for a treatment of 
Lepton Flavour violation in the MSGUT which will be worth exploring in detail if a viable 
region of the parameter space of the perturbative MGUT emerges. On the other hand if no 
such candidate region is available then the remaining options that will need to be explored 
could consist on the one hand of the analogous calculations in the nMSGUT, namely when 
an additional 120-plet is introduced : which radically enlarges the possibilities as far as 
FM fitting is concerned Or the MSGUT could be extended by engineering the models 
to ensure Type II dominance|33- 

Another possibility is that the symmetry breaking at the GUT scale is primarily deter- 
mined not by the perturbative superpotential but rather by a non-perturbative mechanism |3 91 
HO] whereby gaugino condensation in the coset SO(10)/H (where H could be, e.g, SU{5) x 
U{1) or SU{5)) drives a H-singlet Chiral condensate of (say) the 210-plet field at a scale 
~ A/ A" > Mx- In that case the spectra given by us in terms of the vevs a,p,u! are still 
of use but they are no longer determined by the cubic equation^. Rather after breaking 
the symmetry to the group H at a scale higher than the perturbative scale Mx one would 
examine the symmetry breaking in the effective SU (5) symmetric theory at lower energies 
~ Mx- Such a scenario would thus not only utilize the problematic strong coupling regime 
lying just above the perturbative unification scale for dynamical symmetry breaking of the 
GUT theory at a scale Ax ( which is determined by the low energy data and the Grand 
Unified structure and functions as an internally defined upper cutoff for the MSGUT) but 
also provide an explanation for the "SU(5) conspiracy" that seems to operate within 
SO (10) Susy GUTs when proper account is taken of neutrino data, superheavy Susy spec- 
tra and RG evolution. This kind of scenario would fit naturally with a Type I mechanism 
which favours small values of the 126 vev. Thus it contrasts sharply with the proposal 
ofjHZj and should be distinguishable phenomenologically from it. 

We conclude with a tentative proposal to reconcile String theory based models and 
the 126 based Type I and Type II seesaw mechanisms that occur in RGUTs. Recall 
that String theory, particularly with level 1 Kac-Moody algebras, does not favour the 
emergence of effective GUTs containing adjoint and larger representations as its massless 
sector [HI]. Even the use of higher KM levels permits the occurrence of only very restricted 
numbers and types of higher GUT representations. In particular, in the case of SO(IO), 
one cannot enlist the type of combinations characteristic of RGUTs (45 © 54 © 126 © 126 
or 210 ©126© 126) etc. These difficulties led to a decline in attempts to SO(IO) GUTs 
from string theory. However of late the growing appreciation of the naturalness of SO (10) 
unification in the light of discovered neutrino mass [22] has motivated renewed effort in 
this direct ion [HS]. SO(IO) type families are generated with gauge symmetry pre-broken 
to MSSM or somewhat larger. However the issue of implementing the seesaw mechanism 
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whether Type I or Type II in the way achieved so naturally in RGUTs remains problematic 
due to the difficulty of in building models in which the 126 © 126 representations remain 
massless in the string model. 

A way out of this difficulty may perhaps be found by appreciating that in RGUTs the 
126 © 126 fields are superheavy and their neutral components have either a superlarge 
vev (corresponding to Mb-l ~ Mx ) or very small vevs (~ My/ or ~ itli, ). The former 
kind of vev is that of the right handed triplets that give rise to the right handed neutrino's 
superlarge Majorana mass and thus a small Type I seesaw mass for the MSSM neutrino. Its 
large size is compatible with the "pre-broken" structure of "string derived GUTs" where 
the breaking of the SO (10) /GUT gauge symmetry still discernible in the matter super 
multiplet structure is accomphshed already at the level of defining the light sector of the 
String theory and in an effective description corresponds to a vevs of superheavy fields 
in the appropriate large chiral representations (such as 45,54,210,126,126 etc ). On 
the other hand the light vev ~ My/ < Ms > ITeV in MSSM doublet channels occurs 
only after breaking of supersymmetry and its small size is then thought to be naturally 
compatible with the status of Susy breaking as a tiny correction to the pre-broken String 
GUT picture which is supersymmetric, conformal etc and therefore its derivation can be 
legitimately postponed. Furthermore the Type II seesaw mass generating vev of the left 
handed triplet in the 126 is an even higher order effect that arises due to a tadpole in a 
superheavy field induced once EW symmetry breaking has taken place due to the coupling 
of this superheavy field with the doublets that get a small EW weak vev. 

Thus it appears that the 126 © 126 fields need no longer be sought in the massless 
sector of the String theory. Instead it is sufficient to investigate whether the superheavy 
"Higgs Channel" corresponding to 126 or the above mentioned relevant sub-representations 
in fact couple to the putative light fields in an appropriate way. To put it simply the 
implementation of the Rp preserving seesaw mechanism in Stringy GUTs may require a 
small "leakage" connecting the superheavy "field" in the 126 channel to the matter fermions 
and light doublets and the availability of the control parameter Mw/Mx provides a natural 
way to keep the destabilizing effects of such heavy-light couplings under control : given 
that some way has somehow been found to break Susy and generate the EW scale in the 
first place ! This refinement of the effective theory paradigm used to extract the low energy 
theory from string models is both novel and consonant with the characteristic and elegantly 
consistent tying together of very large and very small mass scales achieved by the seesaw 
mechanism. 
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Appendix : Tables of masses and mixings 



In this appendix we collect our results for the chiral fermion/gaugino states, masses and 
mixing matrices for the reader's convenience. Mixing matrix rows are labelled by barred 
irreps and columns by unbarred. Unmixed cases(i)) are given as Table I. 

ii) Chiral Mixed states 



a)[8,l,O](i?i,i?2) = (0/,0;,^i^O)) 



/ (^-^«) -V^^^ (61) 
I -V2Xu} m + X{p-a) I ^ ^ 



rriR^ = l^±l = |2m[l + (| - a) ± ^ + 2u^]\ (62) 
The corresponding eigenvectors can be found by diagonalizing the matrix TZV) . 

b) [1, 2, -l]Chi, h, h, hi) ® [1, 2, l](/ii, /i2, /i3, hi) 

-(Tja yT(15)a ^(15)a 't'l^ \ a, ( fJ V^^^) V^^^) ''^t^S 



n 



—Mh +7\/3(c(; — a) — 7V^(c^ + a) —^cr ^ 

-^^/?,{u + a) -(2M + 4r;(a + a;)) 

7v^(a;-a) -(2M + 477(0 - a;)) -2r]a\/^ 

—(77 — 2?7(T-\/3 —2m + 6A(6i; — ot) y 



The above matrix is to be diagonalized after imposing the fine tuning condition DetTi 
to keep one pair of doublets light. 

c) [3, 1, |] {h: h,h, h: h) e [3, i, ts, is, h) 

= {H^"^, E^*^) , E^^^) , E^o, 04^^+)) e {Hp4^, E(a)^4, E^4(a) , E^4(flo) , <Pfi{R-)) 



Mh 


7(0 + p) lip -a) 


2V2iuj^ 






lip - a) 


2M 










lip + a) 


2M 


4V2icj?7 


2iria 






-AV2iuri 


2M + 2r7j9 + 2r7a 


-2V2ria 




iaj 


2iria 


2^/2ria 


-2m-2X{a+p- 


4u) 1 
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Fie\d[SU{3),SU{2),Y] 


PS Fields 


Mass 


A[l,l,4],A[l,l,-4] 


V2 ' V2 


2(M + 77(p + 3a + 6a;)) 


Ci[8,2,l],C'i[(8,2,-l] 


y A y ^ . 


2(-M + 77(0 + 0;)) 


C72[8,2,1],C2[(8,2,-1] 


A T 

y ■ y ^ 


2(-Af + 77(a-cu)) 


D,\3 2 ^1 DJf3 2 -^1 


E- i E ^• 


2fM + Tjfa + ujW 


D,[3,2,l],D2[{3,2,-l] 


V 4 V 


2(M + r7(a + 3a;)) 


£;i[3,2,|],^i[(3,2,-l] 




-2(M + ?7(a -u;)) 


X[3,1,-|],X[(3,1,|] 


Sp4(ii-),S^R+) 


2(M + 77(a + p + 2u;)) 


L[6,1,|],L[(6,1,-|] 


/■^'(KO) sp'\i,v \ 
MM a/2 


2(M + ?7(p- a)) 


M\6, 1, 11, Mffe, 1, -11 




2(^4" + r)(p - a + 2uj)) 


7V[6,1,-|],7V[(6,1,|] 


/v=r'(-R-) v^'/iP N 


2{M + r]{p - a - 2uj)) 


0[l,3,-2],0[(l,3,+2] 


S44(L) Sf^-j 
\/2 ' \/2 


2{M + r]{3a - p)) 


Pfs 3 _2i pr3 3 21 


^4 


2(M + nia — v)) 


W[6,3,l],W[{6,3,-l] 




2{M -r]{a + p)) 


/[3,l,fl,J[(3,l,-f] 


<^ P(iJ+)' ^4('il-) 


-2{m + X{p + a + 4uj)) 


5[1,3,01 




2{m + X{2a-p)) 


g[8,3,oi 




2{m - \{a+p)) 


C/[3,3,|l,f/[3,3,-|l 




—2{m — X{p — a)) 


F[1,2,-31,V[1,2,3] 


<P44a2 ^al 
V2 ' v/2^__ 

('^JiPai''^aT)M<i> 


2{m + 3X{a + u)) 


E[6,2,|l,5[(6,2,-|l 


-2(m + A(cj - a)) 


F[6,2,-ll,F[(6,2,il 


('^/iPa2'<^aT)M<i> 


2(m - A(a + ou)) 


Z[8, 1, 21, Z[8, 1,-21 




2{m + \{p - a)) 



Table 1: i) Masses of the unmixed states in terms of the superheavy vevs. The SU{2)l 
contraction order is always F'^Fa- The primed fields defined for SU{3)c sextets maintain 
unit norm. The absolute value of the expressions in the column "Mass" is understood. 
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iii) Mixed gauge chiral. 

a)[l, 1, 0](Gi, G,, Gs, G„ G„ G,) ^ {<f>, cl>(^'\ ^"^"y^^"' ) 



/ 



g = 2 



m 





irja 



\ 



V2 





^/^XuJ 

2^/2\uj m + A(p + 2a) 
—i'r]\f?)U 





m + 2\a 







ir]a —zrja 

—iriyf^ia iy/Sija 

M + rj{p + 3a - 6u;) 
M + r]{p + 3a- 6a;) 



b) [3, 2, -|] (^2, ^3, ^4, ^5) © [3, 2, |] (^2, ^3, ^4, ^5) 

= (^4ai' ^('s)a2' ^it^^ Ki) ® (^^2' ^i^l-i ' ^ilvi ' -^Z: 







V5ga* 
2 

2 

J 



^jlaV ^p,4aV ^pAaV 



( -2{M + r]{a - 3uj)) -2^/2i7]a 2ir]a igV2'^ 

2ij2rja -2{m + X{a - u)) -2\f2\u) 2g{a* - lu*) 

-2irja -2V2Xu; -2{m - Xoj) ^/2g{uj* - p*] 

\ -igy/2a* 2g{a* - u*) gV2(u;* - p*) 

c) [l, 1, -2KFi, A, F3) © [1, 1, 2](Fi, F2, F3) 

= (S44(i?o), © (Sj^o)> ^(12)' \r+))- 

/ 2{M + r]{p + 3a)) -2iV3ria -gV2a* 
T = 2i^/3ria 2{m + \{p + 2a)) ^Ugu:*) 

V -gV2a* -\f2iig^* 

d) [3, 1, -|](Ji, J2, J3, J4) © [3, 1, |](Ji, J2, J3, J^) 

= (/)^, (/)/^™\ A/) e (E^4(R+), (Z*^^, (Z'^jIro)' '^a) 



(63) 



(64) 



J 



I 2{M + r]{a + p - 2u;)) 
27]a 
~2V2r]a 
\ -igV2a* 



-2r]a 

-2{m + \a) 
-2^/2\uj 
2^/2iga* 



2V2r]-a 
-2\f2\^ 
-2(m + A(a+p)) 



-ig^o* \ 
-2igV2a* 
-Aigu* 

/ 



e)[3, 2, |](Xi, X2, X3) © [3, 2, -|](Xi, X2, X3) 

-(0S^^0S^^A^t)®('^S.2>£.2>A,,„2) 



(65) 



X 



f 2(m + X(a + oj)) -2V2Xu; -2g(a*+u;*)\ 

-2^/2X10 2{m + Xuj) V2gluj*+p*) 

\ -2g{a* + u;*) V2g{u;*+p*) / 



(66) 
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Figure 1: Plot of the threshold corrections A^*'^\aQ^{Mx)) and A'^^'^\LogiQMx) vs ^ for 
real ^ : real solution for x, A = .7. Note the exclusion of .25 < ^ < 8.6 from the requirement 
that A^^'^^LogioMx) > -1. 
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Figure 2: Plot of the corrections A^^'^^Uq^ (Mx)) and A^^'^\LogioMx) vs ^ for real ^ : real 
solution for x, A = .1. Note how the region around ^ = is excluded by A*^*^)(a^^(Mx)) > 
— 10 while A^^'^^LogioMx) > —1 continues to rule out the 8.6 > ,^ > .25 region. 
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{{NSOL, lam, eta, gam, gambar} , {1, 0.7, 0.5, 0.3, 0.2}} 




Figure 3: Plot of A^^'^^Sin'^OwiMs)) vs ^ for real ^ : real solution for x, A = .7 (upper) 
and A = .01 (lower). Note how the large change in Sin'^9y^{Ms) for 2.8 < ^ < 8 also 
excludes this region. 
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Figure 4: Plot of the corrections A^*'^\aQ^{Mx)) and A'^'^^^LogioMx) vs ^ for real ^ : real 
solution for x, A = .01. Only intermediate values of ^ around the two SU{5) singularities 
are now viable since ac is too large. 
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Figure 5: Plot of the corrections A^^'^\aQ^(Mx)) and A^^'^^LogioMx) vs ^ for real 
complex solution for x, A = .7. Most values of ^ are allowed. 
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Figure 6: Plot of the corrections A^^'^^Uq^ (Mx)) and A^^'^^LogioMx) vs ^ for real ^ : 
complex solution for x, A = .01. Decreasing A further caused a catasrophic increase of 



36 



{{NSOL, lam, eta, gam, gambar} , {2, 0.7, 0.5, 0.3, 0.2}} 



- 



-0.005 - 



-0.01 - 



ti -0.015 - 



^5 -0.02 



-0.025 - 



-0.03 - 



-0.035 -_ 







-0.005 - 



-0.01 - 



-0.015 - 



-0.02 



-0.025 - 



-0.03 - 




-20 -10 10 20 

Xi 



Figure 7: Plot of the threshold and two- loop corrections to Sin'^6w{Ms)) vs ^ for real '■ 
complex solution for x, A = .7 and A = .01. Values of |^| < 10 are allowed. 



37 




Figure 8: Siv?29i2, Sin'^2923 vs ri/r^ for the sample real solution of Bertolini and Malinsky. 
Note the rapid collapse of the large Type II angles as Ita/vs] increases. 
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Figure 9: Plots of MSGUT R = r^/r^ vs ^ for real ^ : real solution(upper) and complex 
solution (lower), A = .7. In the real case the candidate Type II region .8 < ^ < 4 with 
small R is disallowed by A^*'^\LogioMx) > — 1- In the complex case there is no Type II 
candidate region. 
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Figure 10: Plots of the maximum neutrino masses in the MSGUT (obtained using the 
sample solution of Bertolini and Malinsky) vs ^ for real ^ : real solution for x, A = .7. Type 
I (upper) and Type II (lower). The candidate region .25 < ,^ < 8.6 with possibly large 
m^"^ is disallowed by A^^^\LogioMx) > -1. 
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Figure 11: Plots of the maximum neutrino masses in the MSGUT (obtained using the 
sample solution of Bertohni and Mahnsky) vs ^ for real ^ : complex solution for x, A = .7. 
Type I (upper) and Type II (lower). Type I solutions are always dominant and may have 
the right magnitude only in 4 < ^ < 6 region. 
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Figure 12: Magnified plots of tlie maximum neutrino masses in tlie MSGUT (obtained 
using the sample solution of Bertolini and Malinsky) vs ^ for real ^ : complex solution for 
X, A = .7. Type I (upper) and Type II (lower). Type I solutions achieve about 10~^ of the 
required values at their maximum. 
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